Variation of Superconducting Transition Temperature in Hole-Doped Copper-Oxides 
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The experimentally observed difference of superconducting critical temperature T c in hole-doped 
cuprates is studied by using an extended interlayer coupling model for layered d-wave superconduc- 
tors. We show that the change of the maximum T c from series to series is determined by the next 
nearest neighboring hopping t' , while the difference of the maximum T c among the compounds in 
a homogeneous series is controlled by the interlayer pairing strength. Our results provide helpful 
guidelines in the search for new high-T c superconductors. 
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The nature of high temperature superconductors is a 
challenge problem in condensed matter physics. A com- 
mon feature of copper oxide superconductors is the pres- 
ence of CuC>2 plane. It has been observed that the super- 
conducting critical temperature T c varies parabolically 
with the hole concentration nn in Cu02 plane with a 
maximum T™ ax at an optimal doping level Q, Q . How- 
ever, T™ ax attainable is different from series to series, e.g. 
35 K in I^-zSr^CuC^ |3j and 97 K in HgBa 2 Cu04 +( 5 
. An obvious question is what is the crucial parameter 
that governs the T™ ax of each family. 

Among various parameters proposed, the Madclung 
potential at the apical oxygen relative to that at the pla- 
nar oxygens Q was found to correlate with T™ ax rather 
well, pointing to the primary importance of the apical 
oxygens for the electronic structure relevant to supercon- 
ductivity. Further investigations 0, Q revealed that the 
effect of the apical oxygens on high-T c superconductivity 
in reality translates into a correlation between T™ ax and 
the next nearest neighbor hopping parameter t' in the 
t-t'-J model. In these approaches, t' was considered as 
a single parameter reflecting the main difference among 
various cuprates. If we consider the homologous series, 
the universality of such a correlation would be seriously 
questioned. For example, the bilayer and trilayer TI2- 
based and Hg-based compounds have almost same t' [7j, 
but their T™ ax 's are significantly different. 

Our goal in this work is to extract and identify 
which parameters govern the T c behaviors in hole-doped 
cuprates. We apply an interlayer coupling model to CuC-2 
layer systems and then calculate T c based on the BCS 
gap equation with d-wave symmetry. Our results sug- 
gest that the difference of T™ ax from series to series is 
the result of different next nearest neighbor hopping t', 
while the difference of T™ ax between the compounds in a 
homologous series is controlled by the interlayer coupling 
strength Tj. 



The effective layered Hamiltonian we consider is 

Ika Ikk' 
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where Sk is the quasiparticle dispersion, fj, is the chemical 
potential, ct is a quasiparticle creation operator pertain- 
ing to the layer (I) with two-dimensional wave- vector k 
and spin a. The summation over IV runs over the layer 
indices of the unit cell. The intralayer interaction Vkk 1 is 
assumed to be independent of I. The interlayer tunneling 
is parameterized by Tj(k) = Tj(cosk x — cosk y ) 4 Q- 

We assume that the superconducting gap is charac- 
terized by the nonvanishing order parameter b k =< 
c /cT c -fcl > - Based on the BCS theory the gap function 
A l k satisfies the following equation 

Ai = £v fcfc ,&i+r J (A0(&L +1 + &5r 1 ) , (2) 

k' 

where b l k = A l k x l k and the generalized pair susceptibility 
is x{ = (2-Efc) -1 tanh(/3£^/2). The quasiparticle spec- 
trum is E\ = J (e k - m) 2 + |A?c| 2 - 

The spatial dependence of the gap is taken the form 
0: A|. = Afe ±lal . Then the general solution of the 
homogeneous part is 

Ai = A+e^+A^e-^ . (3) 

Considering the fact that the gap vanishes on the layer 
ends / = and n + 1, the natural boundary conditions 
for the gap are A° = A^ +1 = 0. The wave vector of the 
oscillating gap is then determined by 

( e lal e- tal ) ( A k k ) = ° ■ 

The vanishing determinant of the matrix provides a non- 
trivial solution only when a = ^ir/(n + 1) where £ is an 
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FIG. 1: T c vs n H for various t' with V =0.038 eV (a) and 
for various V with t' = —0.02 eV (b) in monolayer cuprates. 



FIG. 2: Calculated T? 



(a) and (b) as a function of t' 



(or J + 2t') for various V in monolayer cuprates. 



integer. Thus we obtain A k = — £ 
of spatial dependence of the gap is then 

£7rm 



At . The solution 



A l k = 2iA k sin 
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(4) 



The solution with the lowest energy is nodeless inside the 
n CuC>2 layers which leads to m = 1 for the supercon- 
ducting state. 

Around critical temperature T c , we can take Xk m 
a simple form: xL — (2i?fc) -1 ta,nh((3Ek/2) = Xk with 



E k = J (e k - n) 2 + |Afc| . Then gap magnitude A k can 
be written as 



A fe - ^ V k k'Xk-A k , = f(n)Tj(k)xkA k 



(5) 



where f(n) — 2 cos(7r/ (n + 1)). 

To account for the experimental observed d-wave gap, 
we assume a ci-wave pairing potential 



Vfefe' = Vg(k)g(k') , g(k) = cosk x - cosk y 



(6) 



The gap magnitude is thus A^ = Aog(k) and the pa- 
rameter Ao is determined by the following self-consistent 
equation: 



1 Vg^k) + f(n)Tj(k) ^ pE k 
M 2- 5 tanh(— ) 



2N 



Ek 



(7) 



The value of T c in layered d-wave superconductors is then 
obtained by solving Eq. (7) at A =0. 

In order to self-consistently calculate T c for a given \x in 
conjunction with the equation determining tih, we need 
an explicit form of e k - It has been established that 
the quasiparticle excitation spectrum of cuprates can be 



well described by the t-t'-J model. Within the framework 
of the t-t'-J model, the dispersion e k is given by fllHl^ 



(J + 2t') cos k x cosk y 



J 



(cos 2k x + cos2k y 



(8) 



For monolayer insulator La2Cu04, experiments [l3j and 
theoretical calculations 0] give a J =0.128 eV. There 
are small variations of J among various Cu-0 insulators 
(lif but we expect a value of J =0.128 eV is a generally 
good representation for all Cu-0 materials. Then one 
can determine T c as a function of njj based on Eqs. (7) 
and (8) once having knowledge of t', V, or/and Tj. 

First we consider the variation of T c in monolayer 
(n = 1) hole-doped cuprates. Figures 1 (a) and (b) 
show the calculated T c in monolayer superconductors as 
a function of uh in some interested parameters range of 
t' and V. As shown, T c initially increases with increas- 
ing nji, takes a maximum around an optimal doping level 
n^* and then decreases with further increasing uh ■ This 
parabolic relation between T c and nu agrees with general 
experimental observations in monolayer cuprates 0, Q] . 
We notice that T™ ax systematically changes with t' , but 
it monotonically increases with V, as one expects. The 
difference between these two parameters is that de- 
pends significantly on t' , while it scarcely changes for 
different values of V. These results indicate that the pa- 
rameters controlling T™ ax would be cither t' or V or both 
of them. 

In Fig. 2 we plotted the t' dependence of both T™ ax 
for monolayer cuprates. As t' increases, T™ ax 
increases and then decreases through a maximum for all 
V studied. The occurrence of the maximum implies that 
the enhancement of T™ ax due to the increase in t' is lim- 
ited. n ^ 1 behaves in a similar manner with t' as TJJ lax . 
For J + 2t' 



and n°fl 



> 0, 



opt 
* H 



decreases with increasing t'. Al- 
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TABLE I: Summary of the experimental results of the critical temperature T™ ax at optimal doping, the distance dc u -0{a) 
between the copper and apical oxygen atoms, the distance dc u ~o{p) between the copper and in-plane oxygen atoms, and the 
calculated values of the bond valence sums of copper Vcu and the difference in the Madelung site potentials AVm for a hole 
between the in-plane oxygen and copper atoms in some typical monolayer cuprates. 

Cuprates T? a * (K) d Cu _ 0(a) (A) d c „-o M (A) Vcu AVm (eV) 

Lai.ssSro.isCuCU 35 2.4124 1.8896 2.539 49.620 

BiaSn.eiLao.sgCuOe+s 36 2.461 1.901 2.437 48.437 

TlBai. 2 Lai. 8 Cu0 5 52 2.500 1.9240 2.280 48.409 

Tl 2 Ba 2 Cu0 6 90 2.714 1.9330 2.135 47.081 

HgBa 2 Cu0 4+a 97 2.780 1.9375 2.091 46.81 



TABLE II: The critical temperature T™ ax and the ratio olTj/V in homogeneous copper-oxides series at optimal doping. The 
brackets are the experimental data taken from the works of Refs. f2ll3ll4ll5l ll6Llr^ .l23ll25ll2^ | . 
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Bi 2 Sr 2 Ca„- 


-lCu TO 2n +4+i 


36.0 (36) 


90.0 (90) 


115.5 (110) 


127.8 


134.7 


150.7 


0.1945 


Tl 2 Ba 2 Ca n - 


_lCu„0 2n +4+<5 


90.0 (90) 


115.0 (115) 


125.2 (125) 


130.1 (116) 


132.9 


139.4 


0.0906 


TlBa 2 Ca„- 


-lCu n 2n +3+J 


52.0 (52) 


107.0 (107) 


131.3 (133.5) 


143.0 (127) 


149.5 


164.6 


0.1930 


HgBa 2 Ca n . 


_lCu„0 2n + 2 +,5 


97.0 (97) 


127.0 (127) 


139.2 (135) 


145.2 (129) 


148.6 (110) 


156.4 


0.1135 



though T™ ax depends on V, is nearly independent 
of V over a wide range of t'. 

To trace the clue to the change of T™ ax among mono- 
layer cuprates^we list in Table I the experimental results 

of T™ ax HHEEII3, the distance d Cu-o(a) between 
the copper and apical oxygen atoms and the distance 
dcu-o(p) between the copper and in-plane oxygen atoms 
taken from the works in Refs. 0,O|, the calculated values 
of bond valence sums (BVS) of copper Vc u and the dif- 
ference in the Medelung site potential for a hole between 
the copper and the in-plane oxygen AVm- To get effec- 
tive BVS of copper, we follow the method proposed by 
Brown 18]. The results of AVm based on the structural 
data are taken from the works in Refs. Here we 

observe one important experimental fact: T™ a:E increases 
systematically with enlarging dc u -o{a)- Band structure 
calculations [Jj revealed that t' increases with d Cu -o(a) 
for the monolayer cuprates reported so far. Thus the in- 
crease of T™ ax with increasing t' should capture the basic 
physics of the monolayer cuprates. 

It has been proposed [2(3, 001 that Vc u and AVm are 
two essential factors governing T c and represent an essen- 
tially equivalent physical content. Materials with larger 
T™ ax tend to have a smaller V Cu or AVm 21]. Since 
the variation of Vn„. or AVm reflects the corresponding 
change of Uh [U |2^, |23| , the increase of the calculated 
rpmax w j^} 1 decreasing n ^ for a wide t' range is obviously 
consistent with the experimental data shown in Table I. 
This n^ 1 dependence of T™ ax is also consistent with the 
muon spin resonance (/iSR) measurements |3]. On the 
other hand, the fact that the change of T™ ax with V 
is almost independent of n ^ 1 (Fig. 2(b)) rules out the 



possibility of V being a dominant factor in governing 
the change in T™ ax . The present results lead us to con- 
clude that the increase of T™ ax with d Cu _ ( a ) among the 
monolayer cuprates is a result of the increase in t' . One 
prediction is that T™ ax decreases with further increas- 
ing t' after through a saturation. Thus, materials with a 
relatively long dc u -o{a) bondlength would not expect to 
have a high T™ ax . 

The values of t' were determined in a self-consistent 
way as follows. From Fig. 2 (a) we learned that there 
exists a maximum for given V . Among the monolayer 
cuprates discovered so far, HgBa2Cu04+a possesses the 
highest T™ ax of 97 K. Assuming this is the highest 
value in all monolayer cuprates, we derived a value of 
V = 0.03762eV from curves of T™ ax versus t'. Equa- 
tion (7) yields t' = —0.0183 eV for the optimally doped 
HgBa2Cu04+a. For other optimally doped monolayer 
compounds with T™ ax < 97K, t' should be smaller than 
—0.0183 eV because of their shorter d Cu -0(a)- The rel- 
ative t' is then obtained by using the experimentally ob- 
served T™ ax . 

Next we consider n, the number of Cu02 layers, de- 
pendence of T c in the layered homogeneous series. In 
general, T™ ax initially increases with n, maximizes at 
n = 3, and then decreases with further increasing n |25l |. 
To calculate T c for multilayers, we use the same disper- 
sion £fc and V as obtained from the monolayer. The 
interlayer tunneling strength Tj is determined by using 
the experimental values of T™ ax for monolayer and bi- 
layer compounds in the same homogeneous series. As an 
example, in Fig. 3, we show curves of calculated T c ver- 
sus riH as a function of layer number n in the Hg-based 
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ation in hole-doped cuprates on the basis of an extended 
interlayer coupling model. We demonstrate that the next 
nearest neighboring hopping t' dominates the variation of 
the maximum T c from series to series and the interlayer 
coupling strength controls the difference of the maximum 
T c among the compounds in a layered homogeneous se- 
ries. These results provide helpful guidelines in the search 
for new high-T c superconductors. 

The authors are grateful to J. S. Schilling and W. G. 
Yin for many helpful discussions. This work was sup- 
ported in part by the Earmarked Grant for Research of 
Project CUHK 4037/02P. 



FIG. 3: Calculated T c vs nn in HgBa2Ca n _iCu„02n+2+* as 
a function of the number of G1O2 layers. 

series. The parabolic behavior is generally observed for 
any layered compound. The calculated T" lax in four typ- 
ical homogeneous series are summarized in Table II. The 
experimental results are also listed for comparison. As 
can be seen, T™ ax initially increases with increasing n 
and then saturates as n — > oo. This behavior is in good 
agreement with those obtained from both the interlayer 
mechanism (j| H(| and Ginzburg-Landau theory j^, • 
The upper limit of T™ ax for infinite layer compound is in 
the range of 139.4 to 164.6 K. The highest T™ ax of 164.6 
K is found in the Tl-based series. Our results for n = 3 
agree with experiments very well. The predictions made 
here for T™ ax of the trilayer compound is the best ones 
compared to previous theories [Si |26j, |2jj, |2£| • 

The present study shows that interlayer coupling is 
the driving force for the enhancement of T™ ax for mul- 
tilayer systems. This does not conflict with the exper- 
imental fact that T™ aa: saturates as n > 3. In fact, 
there exist five-fold (outer) and four-fold Cu02 (inner) 
planes surrounded by pyramidal and square oxygens in 
the multilayer system. Investigations carried out by dif- 
ferent experimental techniques and model calculations 
[23L 12^. lifl l3l| showed that the distribution of charge 
carriers are nonhomogeneous among the Cu02 sheets and 
the hole concentration in the outer Cu02 plane is larger 
than that in the inner Cu02 plane. BVS analyses |23j 
and NMR studies :3lj on the Hg-based series revealed 
that the highest T™ ax corresponds to the smallest differ- 
ence in tih between two types of Cu02 planes. When the 
number of Cu02 layer is larger than three, the reduction 
of T™ ax comes from the large difference in njj between 
the outer and inner Cu02 planes. For compounds with 
more than three Cu02 planes, the enhancement of T™ ax 
seems possible at ambient pressure if one can adequately 
dope the inner planes. 

In summary, we have investigated the observed T c vari- 
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